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Using the EAS size spectra measured with the MAKET ANI array on Mt. Aragats, 
Armenia (3200m a.s.L- 700g-cm~ 2 ) in the range of N e = 10 5 — 10 7 for different angles- 
of-incidence, the EAS attenuation length has been determined applying different 
analysis methods. The analysis is based on a data sample of 2.5 • 10 6 events collected 
in the period of June, 97 - April, 99. The results are compared with results deduced 
from data of the EAS TOP and KASCADE experiments. 
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1 Introduction 



The intensity of Extensive Air Showers (EAS) with fixed shower sizes N e is assumed to decrease 
exponentially with increasing atmospheric depth of the observation level. This is considered to be 
due to the absorption of the particles of the EAS cascade following an exponential law 
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N e (X) = N e (X )expi^ j-^-j, withA>AV (1) 
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An is a definite initial atmospheric depth after the maximum of the longitudinal development where 
the number of (charged) particles is N e (Xo) and further decreasing exponentially, N e (X) is the num- 
ber of particles of the EAS at the slant depth X[g ■ cm~ 2 }. 

The quantity A controls the attenuation of particles of the individual cascade [J (size attenuation 
length). It is related to the inelastic cross sections (to the mean free path length A a) of the interaction 
of the primary cosmic ray particles with air nuclei. The attenuation of the flux intensity of Exten- 
sive Air Showers is characterized by a related quantity Xn (intensity attenuation length, absorption), 
which can be directly measured by cosmic rays detector arrays. Thus measurements of the attenua- 
tion of the EAS intensity in the atmosphere are considered to be an interesting source of information 
about hadronic interactions, especially if extended to the ultrahigh energy region expected from the 
forthcoming LHC and TESLA accelerators. In addition due to the sensitivity of the cross sections 
to the mass of the primary, alterations of the attenuation length with the energy may be indicative 
for the variations of the mass composition. Measured results imply tests of the energy dependence 
of the extrapolated cross sections used for Monte Carlo simulations. 

The investigations of the present paper are based on an EAS sample measured 1997-1999 with the 
MAKET ANI array [21 Ej on Mt. Aragats station (Armenia) and registered for different angles-of- 
incidence in the zenith angle interval = — 45° . The data basis of the analysis can be enlarged by 
published data from KASCADE (1046 g-cm~ 2 ) 4 and EAS TOP (810 g-cmT 2 ) experiments. Spec- 
tra measured by EAS TOP are given in Ref. 6 . Data and zenith angle dependence for KASCADE 
results are obtained by scanning the spectrum plots communicated by the KASCADE collaboration 

m 

We apply different procedures to deduce the attenuation. First we consider the degradation of the 
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EAS flux with fixed shower size N e with increasing zenith angle i.e. increasing atmospheric thickness 
of the shower development (characterized by the intensity attenuation length (Ajv) |Hj)- Differently 
the technique of the constant intensity cut (CIC) t 9„ considers the intensity spectrum of EAS events 
and relates equal intensities observed at different atmospheric depths. 

There is the tacit assumption that the shower size reflects the energy of the primary. The procedure 
can be refined by using the knee position in the N e spectrum as a bench mark for a well defined energy, 
so far we may associate the knee phenomenon to a feature of the primary energy spectrum of cosmic 
radiation. 

2 Experimental spectra 

The experimental basis of the present investigations are measurements of shower size spectra in the 
knee region and their zenith-angle dependence performed with the MAKET ANI array of the Mt. 
Aragats Cosmic Ray Station (3200 m a.s.l.) in Armenia. Details of the measurements and the ex- 
perimental procedures taking into account the detector response are given elsewhere |10] 111]. For a 
detailed description of the knee region the traditional approximation with two different spectral in- 
dices below and above the knee, defining the knee position as intersection of two lines in a logarithmic 
presentation, appears to be insufficient. Hence a more sophisticated method has been applied with 
parameterization of the slope of the spectra (see Ref . [12] ) . 

Tab.l compiles the characteristics of the size spectra measured with the MAKET ANI installation, 
the changes of the slopes in the knee region (AiV efc ), expressed by different spectral indices below (71) 
and above (72) the knee position N ek for the zenith-angle range of = — 45°. For the display and 
the analysis of the zenith-angle dependence, the size spectra are determined in 5 angular bins of equal 
AsecB widths. The accuracy of the zenith angle determination is estimated to be about 1.5° |l()j . 
A correction due to barometric pressure changes, which lead to small fluctuations of the atmospheric 
absorption, has not been made. Figure ^ displays the spectra of mean values of each atmospheric 
depth bin and compares with the results from EAS-TOP 6 and KASCADE |7] experiments. 



I(10 5 <N e < 1.15 • 10 6 ) (8.95 ± 0.18) • 10~ 11 (iV e /10 5 )' ri 



I(N e > 2.56 • 10 6 ) (3.23 ± 0.40) • 10 -13 (iV e /10 6 )'> 2 

71 -2.54 ±0.012 

72 -2.94 ± 0.042 

A(iV e J (1.15 ± 0.034) • 10 6 - (2.56 ± 0.063) • 10 6 

N ek (1.75 ± 0.05) • 10 6 

I(N ek ) (5.83 ± 0.14) • 10~ 14 



Table 1: Flux [m 2 s 1 sr x ] and knee region parameters of the size spectra measured with the MAKET 
ANI array. 

Following fixed intensities of the experimental spectra (see sect. 3. 2) the average N e cascade devel- 
opment can be immediately reconstructed as shown in Figure |2j Note that the results in the range of 
the slant depth observed with the ANI array deviate from the exponential decrease (eqU). That is an 
interesting feature which can be revealed more clearly when combining spectra accurately measured 
on different altitudes. In the present paper we base the formulation of the procedures estimating 
the attenuation on the exponential decrease (eqU>. It is our interest to explore, if this assumption 
applied to the ANI and KASCADE data lead to consistent results. 




Figure 1: Differential size spectra for differ- 
ent zenith angles ranges observed with MAKET 
AN I array, compared with spectra reported by the 
KASCADE Q/ and the EAS TOP [6] collabora- 
tions. 



Figure 2: N e cascade in the observed range of 
the atmospheric slant depth. 



3 Procedures for inference of the attenuation length 
from size spectra 

We consider the differential and integral size spectra I(N e ,X) and /(> N e ,X), respectively. In 
addition to the basic assumption of exponential attenuation of N e (eqjj) a power-law dependence of 
the size spectrum 

I(N e ,X)(xN e ~\ (2) 

with the spectral index 7 is adopted. 



3.1 Attenuation of the intensity of fixed iV e : absorption length 

For different fixed values of shower size N e , on different depths in the atmosphere or/and different 
zenith angles of incidence, from measured spectra (see vertical dotted lines on Figure OJ) we obtain 
several values of corresponding intensities from the equivalent depths from 700 till 1280 g ■ cm -2 . 
Fitting the depth dependence of the intensities by the straight line (in logarithmic scale) according 
to equation: 

I(N e ,X) = I(N e ,X )exp (3) 

we obtain the estimate of the absorption length A at. The absorption length can be estimated both 
by integral and differential spectra. 




Figure 3: Integral size spectra for different zenith 
angles ranges observed with MAKET ANI array, 
compared with spectra reported by the KASCADE 
illustration of the procedures for absorption 
and attenuation length estimates. 



Figure 4: The variation of the knee position with 
the atmospheric depth. 



3.2 Constant intensity cut 

The basic idea of this procedure is to compare the average size of showers which have the same rate 
(showers per m 2 ■ s ■ sr) in the different bins of the zenith angle of shower incidence and different slant 
depth, respectively [§]. 

Considering two different depths in atmosphere Xi,X 2 > Xq the expressions of differential intensities 
I(N e ,X) has the form 
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With simple transformations we obtain: 
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(5) 



The attenuation lengths, obtained by integral spectra do not depend explicitly on spectral index: 

X 2 — ^1 



Aint(I) 
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(6) 



Practically the estimate of the attenuation length is obtained by fitting the N e dependence on the 
depth in atmosphere by the straight line according to the equation The sequence of N e values 
is obtained according to the fixed values of the flux intensity, selected from the interpolation of the 
differential or integral size spectra. 



For each N e value, the slope index 7 used in equation [SJ is obtained by averaging over all used slant 
depths. Selecting equal intensities (~ primary energies) corresponding to different shower sizes N e 
and different depths the value of A^j-j(I) is estimated. Intensity values from 10 -9 to 5. • 10 -6 were 
used for CIC method. 



3.3 Attenuation of the size of the knee 

A special variant of the constant intensity cut is to follow the decrease of the shower size at a constant 
primary energy in the size spectrum. Assuming that the knee phenomenon reflects a feature of the 
primary flux, the variation shower size at the knee with the zenith angle provides the possibility to 
extract the attenuation length. 

Considering the assigned knee position of the data from various experiments, differences within 30% 
are noticed for all X-bins. 

The knee positions obtained by the differential and integral spectra are a bit shifted to the smaller 
N e values (see Figure HJ). The shift is approximately uniform over all investigated depths interval, 
therefore the estimates of the attenuation length by the differential and integral size spectra are very 
close to each other. 



3.4 The relation between the absorption and attenuation length 

We consider the quantity I(N e , X)dN e - the number of EAS at the depth X which comprise N e to 
N e + dN e particles: 

I(N e ,X)dN e ~ N^exp - (7 - 1) X ~ A X ° dN e (7) 

With eqOlwe obtain: 

Knff(N e ) = ( 7 (Ay - 1)X N , (8) 

where, j(N e ) is the differential size spectra index (here we indicate the N e dependence of the slope 
index explicitly). For the integral spectra: 

^int(N e ) = ~f(N e )\ N , (9) 

where, 7(A r e ) is integral size spectra index. 

For the evaluation of the inelastic cross section and for comparison of the three methods described 
above we propose to use the calculated values of the attenuation length A (instead of using absorption 
length Ajv). The attenuation of the number of particles in the individual cascade is more directly 
connected with the characteristics of the strong interaction and is independent from the parameters 
of the cosmic ray flux incident on the atmosphere. In turn the absorption length, i.e. the attenuation 
of the CR flux intensity, reflects also characteristics of the primary flux and is dependent on the 
change of the slope of the spectra. 



3.5 Estimate of the inelastic cross section 

The inelastic cross sections, of the primary nuclei with atmosphere nuclei is related by |Hj: 

2 41 • 10 4 

aTUmbarn) = ^ ^ (10) 



where A denotes the primary nuclei. The quantity is the interaction length of the A-nucleus 
in the atmosphere (note: in some publications the interaction length is denoted by Ajv, where N is 



primary nuclei, in contrast in this paper N is reserved for the shower size). The interaction length 
\a is related with the absorption length by 

\ A = K{E) ■ A A (11) 

The K(E) coefficient reflects peculiarities of the strong interaction model used for simulation. The 
value of the parameter K has to be determined by simulations of the EAS development in the 
atmosphere. Such studies require the development of procedures for the selection of EAS initiated 
by primaries of a definite type (see for example in |131 114] ^ . 

4 Application to the data 

The mean values of the attenuation lengths obtained by various methods from data of the ANI and 
KASCADE installations, as well as for the joint ANI & KASCADE data by the differential (Adi//) 
and integral spectra (A m t) are compiled in the Tables 12^3141 



Min. depth 


MAKET ANI 


ANI+KASCADE 


KASCADE 


X ,g ■ cm' 2 


Ami 


Adi/ 


A m t 




A m i 




700 


248 ± 27 


247 ± 42 


203 ± 10 


203 ± 13 






758 


236 ± 32 


237 ± 51 


195 ±8 


196 ± 12 






816 


211 ±43 


218 ± 70 


186 ±9 


188 ± 13 






1020 










181 ± 14 


182 ± 23 



Table 2: Attenuation lengths for the data from the MAKET ANI and KASCADE installations esti- 
mated by the CIC method from differential and integral size spectra 



Min.depth 


MAKET ANI 


ANI±KASCADE 


KASCADE 


X ,g ■ cm' 2 


A m i 


Arfj/ 


Ami 


Adi/ 


Ami 


Adi/ 


700 


239 ± 14 


240 ± 15 


191 ± 11 


193 ± 13 






758 


232 ± 13 


228 ± 19 


186 ± 10 


184 ± 17 






816 


213 ± 14 


219 ±27 


179 ± 11 


181 ± 24 






1020 










181 ± 7 


183 ±11 



Table 3: Attenuation lengths for the data from the MAKET ANI and KASCADE installations esti- 
mated by the recalculation from the absorption length for differential and integral size spectra 



Min.depth 


MAKET ANI 


ANI±KASCADE 


KASCADE 


X ,g ■ cm- 2 


Ami 




Ami 


Adi/ 


Ami 


A^/ 


700 


302 ± 71 


295 ± 83 


241 ± 17 


237 ± 15 






758 


272 ± 51 


263 ± 42 


242 ± 20 


221 ± 17 






816 






225 ± 21 


225 ± 19 






1020 










232 ± 26 


222 ± 28 



Table 4: Attenuation lengths for the data from the MAKET ANI and KASCADE installations, 
estimated by the "attenuation of knee position" method from differential and integral size spectra 

The alternative estimates of the attenuation length reflect the inherent uncertainties of the methods 
and the statistical errors, as well as the fluctuations of cascade development in the atmosphere, the 
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Figure 5: Attenuation Length dependence on 
Spectra Intensity (Primary Energy). 



Figure 6: Attenuation Length dependence on the 
Shower Size N e . 
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energy dependence of the inelastic cross sec- 
tion and possible changes in mass composi- 
tion. As obvious in Figure[2 the values corre- 
sponding to the minimal equivalent depths of 
used MAKET ANI data, deviate significantly 
from the exponential decrease. The obser- 
vations reflects the flattening of the cascade 
curve just after the shower maximum in the 
altitude 500 — 600 g ■ cm~ 2 . Due to these fea- 
tures the attenuation lengths calculated by 
MAKET ANI data appear to be significantly 
larger than those derived for the KASCADE 
data (Tables El) • 

Therefore, for the combined analysis of the 
KASCADE and ANI data we omitted the 
first and the second zenith angle bins of 
MAKET ANI and calculate the attenuation 
lengths by the remaining 9 (minimal equiva- 
lent depth 758 g-cm~ 2 ) and 8 (minimal equiv- 
alent depth 816 g ■ cm~ 2 ) angular bins. The 
dependences of estimated values of attenua- 
tion length on the shower size and flux inten- 
sity for different amount of the angular bins 
used, are displayed in Figures (note, that 

higher intensities on the X axes correspond to the lower primary energies) and El The attenuation 
length estimates obtained from the differential and integral spectra agree fairly well. The results of 
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Figure 7: Attenuation length obtained by joint analy- 
sis of the MAKET ANI and KASCADE data. 



both CIC and recalculation from absorption length agree within the error bars. The results obtained 
by the "attenuation of knee position" are larger for MAKET ANI and KASCADE. As pointed out by 
S. Ostapchenko ^S] it is the consequence of the large EAS fluctuations with the tendency to shift the 
knee position to the lower energies (and correspondingly to higher fluxes) in a way to "slow down" 
the cascade curve attenuation. 

Well below the shower development maximum starting from 816g • cm -2 KASCADE and MAKET 
ANI data could be fitted with one decay parameter (see Figure EJ) . There is a concentration of the 
knee positions on the curve showing the dependence of the attenuation of the flux intensity (~ pri- 
mary energy). In turn, the curve displaying the dependence of the attenuation length on the shower 
size demonstrates a rather large dispersion of the "knee positions". These observations in size and 
energy scales may be interpreted as an indication of the astrophysical nature of the knee phenomenon. 

5 Conclusion 

Experimental studies of EAS characteristic like the depth of the shower maximum X max , the elonga- 
tion rate dX max / dlogi^E and the attenuation length A are of particular importance, since they map 
rather directly basic features of the hadronic interaction. Strictly, however, the interpretation of these 
quantities in terms of hadronic cross sections cannot bypass the necessity of detailed calculations of 
the shower development. Nevertheless these type of EAS quantities, if compared with Monte Carlo 
simulation results, provide stringent tests of the interaction model ingredients of the simulations. 
The recent results of various experimental installations are sufficiently accurate to enable relevant 
studies of this kind, and combining the data from arrays situated on different altitudes (like MAKET 
ANI and KASCADE) allows a large span in the atmospheric slant depth for reconstructing the de- 
velopment of the charge particle size. In fact such studies, if using a sufficiently large data sample, 
could be continued in a more detailed manner by separating the muon component and taking into 
account the deviations from the exponential shape of the cascade decline. The penetrating muon 
component contributes with smaller attenuation to the development of the considered charged par- 
ticle component, but hardly with an exponential degrading (according to eqQ). Actually by use of 
methods in progress to isolate different primary groups ("pure nuclear beams") of the size spectra 
^UE3> these kind of interaction studies would get of extreme interest. 
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